Two-fluid description of magnetic excitations in iron pnictides 
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We present a phenomenological, two-fluid approach to understanding the magnetic excitations in 
Fe pnictides, in which a paramagnetic fluid with gapless, incoherent particle-hole excitations coexists 
with an antiferromagnetic fluid with gapped, coherent spin wave excitations. We show that this two- 
fluid phenomenology provides an excellent description of NMR data for magnetic "122" pnictides, 
and argue that it finds a natural justification in LSDA and spin density wave calculations. We 
further use this phenomenology to estimate the maximum renormalization of the ordered moment 
that can follow from low-energy spin fluctuations in Fe pnictides. We find that this is too small to 
account for the discrepancy between ab intio calculations and neutron scattering measurements. 
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The discovery that, suitably doped, Fe pnictides can 
superconduct at temperatures greater than 50K— has 
sparked a sudden rush of interest in these materials. 
Much of this attention has been focused on the mag- 
netism of the parent compounds, which is widely believed 
to play an essential role in their superconductivity 2 - - — . 
To date, most theoretical approaches to this problem 
have stressed either the itinerant nature of electrons in 
Fe pnictides^—, or used strong electronic correlation to 
justify mapping them onto a frustrated local moment 
models**. In this paper we embrace the fact that Fe 
pnictides are both metals and magnets, proposing a sim- 
ple, phenomenological, two- fluid description of their mag- 
netic excitations. We argue that spin excitations at low 
energies and temperatures are dominated by the gap- 
less, incoherent particle-hole excitations, characteristic of 
a metallic paramagnet, while for energies and tempera- 
tures comparable with a spin gap A CT , coherent, collec- 
tive excitations of the magnetic order come into play. It 
follows naturally from experimental and theoretical de- 
terminations of the band structure that these two fluids 
are essentially independent. 

While this two-fluid phenomenology is not dependent 
on the underlying microscopic model, it finds a natu- 
ral justification in recent spin density wave calculations, 
which assign metallic and magnetic electrons to differ- 
ent sheets of the Fermi surface^— . It is also consis- 
tent with the band structure measured by photocmission 
studies and calculated by LSDA^. Here we show that 
it provides an excellent description of NMR experiments 
on Fe pnictides with 122 structur o 16 i 17 . We also criti- 
cally re-examine attempts to understand magnetism in Fe 
pnictides in terms of frustrated local-moment models^—. 
While a two-fluid phenomenolgy does not rule out frus- 
tration per se, we find that quantum fluctuations cannot 
account for the observed reduction of the ordered mo- 
ment relative to LDA calculations^. 

Both the magnetic and metallic properties of Fe pnic- 
tides originate in outer-shell Fe 3d-electrons. Band struc- 
ture calculations^*^, supported by photoemissio n 19 ! 20 
and quantum oscillatio n 21 ! 22 experiments, suggest that 
these hybridize with As 4p orbitals to form a Fermi sur- 



face with two electron-like and three hole-like pockets, 
when viewed in a "natural" unfolded Brillouin zone based 
on Fe sites. The fact that magnetic Fe pnictides are met- 
als implies that some part of this complex Fermi surface 
remains gapless, and will support incoherent particle-hole 
excitations with vanishing energy. We treat this as the 
first of our fluids, characterized simply by an average 
density of states at the Fermi energy, no, which can be 
estimated from heat capacity measurements. 

Neutron scattering experiments 2 ^—, meanwhile, re- 
veal a commensurate, collinear, antiferromagnetic (AF) 
ground state with ordering vector k* = (ir, 0,7r), and 
ordered Fe moment w l/^s, much smaller than pre- 
dicted by ab initio calculations^. A single branch of 
spin wave excitations with dispersion, 



LU k , = VA 2 + (v.k') 



(1) 



is found above a gap A CT w lOmeV at the ordering 
vector k' = k — k* = (0,0,0). Spin wave velocities 



{v x ,v. 



are anisotropic, with v x 



> v z . The 



collective excitations of this magnetic order form our sec- 
ond fluid, and, following [27j , we characterize them using 
a quantum non-linear sigma model, 



S = 



1 



2abc 



dxdt [H 2 X ±(d t n) 2 - p x (d x ii) 2 - PyidyXif 



-p z (d z iif + x±Kn x ] , 



(2) 



where \± is the static perpendicular susceptibility, p x , p y 
and p z are spin stiffness' along the major crystal axes 
a, 6, c, and A 2 is an easy axis anisotropy. 

For A (j — > 0, this action describes the long- wavelength 
Goldstone modes, which follow from the symmetry of 
the magnetic order. As such, it can be derived from 
any microscopic model that respects these symmetries, 
whether localized or itinerant. For finite anisotropy 
A CT > 0, Eq. @ predicts a gapped, two-fold degener- 
ate cone of spin wave excitations with exactly the form 
Eq. ([1]), where v a = \J p a /x±- Within a spin density 
wave picture, Eq. ([2]) should remain valid up to an en- 
ergy scale of the spin-density wave gap, estimated to be 
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TABLE I: Leading temperature correction to the ordered moment 
(m s = mo — 8m$ ~ & m s{T)) from spin wave excitations. The 
form of corrections depends only on the gap, A CT . The prefactor 
is determined by the geometric mean spin wave velocities, = 
v x VyV z [(meVA) 3 ], bare moment mo [a*b], transverse susceptibility 
Xx [meV -1 ] and Fe-Fe lattice parameters a,b,c [A]. 



Asdw ~ 31meV for LaFcAsGi^. For the specific case of 
CaFe2As2, it breaks down at energies of approximately 
150meV, where the spin wave branch is seen to enter a 
continuum of excitations 2 ^. 

Our final approximation is to ignore all coupling be- 
tween these two fluids. This appears justified for two 
reasons. Firstly neutron scattering studies^— observe 
sharp cones of low energy magnetic excitations with no 
evidence of the damping that would be expected if the 
spin waves could scatter from the metallic fluid. Sec- 
ondly LSDA calculations^, photoemission studies! and 
SDW theory-^ 3 -"— show no evidence, in the ordered state, 
for a nested pair of Fermi surfaces with spanning vector 
(7r, 0). Hence, within a band picture, there are no avail- 
able particle-hole states close to the ordering vector for 
the spin waves to decay into. 

A simple test of our phenomenology is provided by the 
temperature dependence of the ordered moment Sms(T). 
Within the two-fluid picture, this is controlled by the 
thermal excitation of spin waves, as described by Eq. (|2|). 
Table U summarises our predictions. At temperatures 
relevant to experiment, the spin gap dominates, and in 
Fig. [1] we compare the predicted form of 5ms{T) with 
the low temperature ordered moment, as measured by 
NMR experiments on BaFe2As^ and SrFe2As2^1. We 
have checked that similar fits can be made for 5ms(T) 
obtained from [iSR for LaOFeAs 2 ^ and SrFcaAs^. 

NMR experiments also probe spin excitations through 
the nuclear spin lattice relaxation rate, 1/Ti. This 
has been measured for As nuclei in BaFe2As2^ and 
SrFe2As2^£. For hyperfine interactions, the relaxation 
rate is given b y 31 ' 32 , 



q 



2 x"(u ,q) 

Hun 



(3) 



where is the gyromagnctic ratio of the nucleus in 
question and ujq is the nuclear excitation energy. Both 
fluids contribute to 1/Ti, but at low temperatures the 
leading contribution will come from gapless particle-hole 
pairs within the paramagnetic fluid. We assume that 
the hyperfine interaction is roughly constant, |^4 q | 2 ~ 
|.4o| 2 [(T '/ '/xb) 2 ], over the relevant sheet of the Fermi sur- 
face, in which case the contribution to 1/Ti will be linear 
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FIG. 1: (Color online). Temperature dependence of the or- 
dered moment <5mg(T) as determined by NMR measurements on 
BaFe2As2^ (blue circles) and SrFe2As2i£ (black squares). Data 
is plotted as log \8ms(T)/T^ J vs A^/ksT, where the values of 

Asa = H4Jf and A$ r = 75K are taken from inelastic neutron 
scattering experiments2iL2&. Straight lines show the expected form 
of corrections at low temperatures. The intercept gives the prefac- 
1.6 x 10~ 4 y, B K~§ and M^-f a 5 x 1CT 5 ^i B K~i 



tors M^ a 



in T m 



1/Ti nc - « -nh^A^gl (abcfn 2 k B T = C inc T, (4) 

where gh is the Lande g-factor and we use units where 
A's = 1. 




FIG. 2: (Color online). Simultaneous fits to nuclear relaxation 
rate Tj~ 1 and sublatticc magnctsation mg(T) for BaFe2As2^ 
(blue dots) and SrFe2As2i£ (black squares). The dashed lines 
show the contribution of incoherent particle- hole pairs Eq. (2); the 
full lines show the combined fit including the contribution of co- 
herent, thermally-activated spin waves Eq. 10. The gap values 
Ago = 114-ff and A$ r = 75K are taken from inelastic neutron 
scattering experiments2£i2ii. Insets show simultaneous fits to the 
sublatticc magnetisation, mj. The prefactors are determined to 
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At higher temperatures, the Raman scattering of ther- 
mally excited spin waves also plays a role in nuclear spin 
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Fit 


1.6 x icr 4 


0.032 


7.5 


Estimate 


(0.15 - 2.2) x 1CT 4 


~ 0.015 


0.14 - 33 



TABLE II: Quantitative analysis of 5m s (T) and T, 1 in BaFe2As2. 
We determine the prefactors M^ a (Table Q}, Cf£ c (Eq. Q, and 
C coh ( E 1- EJ by fitting NMR experiments (cf. Figs. [Q [2]l and 
compare these with the quantitative estimates which follow from 
known values of the hypcrfine interactions Aq ~ 1.88 T/fiB and 
Bac = 0.43 T/fiB^, spin-wave velocities 95 < v s < 228 meV X2^, 
spin-gap A CT = 9.8(4) meV — , ordered moment mo = 0.87 Hb—^i 
perpendicular susceptibility Xi = 10 — 4 emu/mol—, lattice con- 
stants [a.b,c] = [2.80,2.79,6.47] Ai£, and density of states n = 
5.8 X 10 24 J-im-Sia. 



relaxation. NMR probes the longitudinal susceptibility 
x'u (u>o, q), which can be calculated directly from Eq. ((2|). 
We consider relaxation due to hyperfine coupling of the 
As nucleus to a four site plaquette of nearest neighbour 
Fe sites, in the same spirit as (33[. For external fields 
applied in the FeAs plane we calculate the form factor, 



frustrated exchange interactions, where the ordered mo- 
ment is strongly renormalized by quantum fluctuations^. 
By analogy, it has been suggested that magnetic excita- 
tions in Fe pnictides can also be understood in terms of 
a frustrated local-moment model, 

"H = J\ x Sj.Sj + J\ v Sj.Sj + J\ z Si.Sj 

(ij)lx (i3)ly (ij)lz 

+J 2 £ SoSj K xy m 2 (Sf) 2 ) + K z £(S*) 2 (7) 

<ij>2 i » 

where {ij)\ a counts first-neighbor bonds in the a- 
direction, (ij)2 second- neighbour bonds in the x-y plane, 
and K xy is a single-ion anisotropy. It is interesting, there- 
fore, to ask what constraints our two-fluid phenomenol- 
ogy places on this effective local-moment picture ? 

A telling, and direct, comparison can be made in the 
context of the ordered moment. At a mean field level, 
the collinear "stripe" phase of the square-lattice J1-J2 
Heisenberg model becomes unstable for J% <\J\\/2 37 1. 



I A, 



42? 2 C (1 - cos q x + cos q y - cos q x cos q y ) , (5) 



where B ac = 0.43 T/fis [3 is the hyperfine matrix ele- 
ment relevant for Raman relaxation processes. The imag- 
inary part of the susceptibility is very sharply peaked for 
scattering by the ordering vector, and hence we approx- 
imate the slowly varying form factor by |-4 q | 2 = 16S^ C , 
and sum over a small region surrounding q = (tt,0,tt). 
We find, 
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coh. 
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(6) 



and 



where = x 2 hi 1 (e~ 1 ^ x ) + x 3 hi 2 (e~ 1 / x ) 1 

Li n (z) = X)fe°=o z k /k n is the n th polylogarithm of z. 

We are now in a position to compare directly with ex- 
periment, and, in Fig. [2j we show the results of simultane- 
ous fits to NMR data for 8m s {T) and 1/T X in BaFe 2 As2^ 
and SrFe2As2^. We treat the total relaxation rate as 
the sum of the contributions of the two fluids, Eq. (j4| 
and Eq. © and fit the prefactors M^ T a (Table [J), Cj^ 
(Eq. 21), and C^ h (Eq. [6]), using the experimental value 
Asa = HAK— for the gap. The agreement is excellent. 
It is also possible to make independent, quantitative es- 
timates of these fit parameters by substituting known 
values of the model parameters directly into the relevant 
equations. In Table [TT1 we show that these two approaches 
are quantitatively consistent for BaFe2As2^ 

One of the important issues in Fe pnictidc magnetism 
has been the size of the ordered moment mg. Fe and 
its oxides typically show a large ordered moment at low 
temperatures. First principles calculations for magnetic 
Fe pnictides suggest that tos«1.5-1.7^^£. The mo- 
ment measured by neutron scattering, in contrast, ranges 
from 0.25hb (NdFcAsO)^ to \\x B (SrFe 2 As 2 )^. The 
AF "stripe" order found in Fe pnictides has also been 
observed in quasi-two dimensional insulating oxides with 



Approaching this transition, quantum corrections to the 
ordered moment diverge, as illustrated in Fig. [3J For AF 
Ji, the dominant correction to mg comes from spin waves 
near the ordering vector. These are described by Eq. 
with v, = A a = 0, and we find, 



Smg = 



mo 



dk 



2x_l (2tt) 2 7| k | <A w k 4tt 2 x_lu 3 



Ki(k), (8) 



where A is a momentum cut-off reflecting the size of 
the spin- wave "cone", Ki is a complete elliptic inte- 
gral of the first kind, and k = y/l — (v y /v x ) 2 . At the 
limit of the (n, 0) AF phase, v y — > 0, and 5mg diverges 
logarithmically^. The contribution of spin waves at 
higher energies must be determined separately, but for 
present purposes can be approximated by a constant off- 
set PS 0.1/!_B. 

At first sight, fine-tuning a J1-J2 model into a re- 
gion with v y <C v x offers the possibility of achieving 
any desired renormalization of the ordered moment, mg, 
cf4~— . The same would hold of any itinerant electron 
model which could be mapped onto Eq. ([2]). However, 
neutron scattering results for Fe pnictides suggest that 
v y ps v x pBI |. Moreover, they clearly show a spin gap 
A CT , and out-of-planc dispersion v z , both of which act to 
cut-off the divergence in 5mg. 

For a gapped, three-dimensional dispersion with 
v y > y/v x v z (A/ir) 3 , Eq. © predicts 
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(9) 



For parameters relevant to BaFe2As 2 , with 
A ps 0.2n/a [2(|, this implies Smg ps 0.13/is, a value too 
small to explain the gulf between ab intio calculations 
and experiment, although comparable with the smaller 
discrepancy with model based SDW theory^. 
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FIG. 3: (Color online). Zero temperature sublattice magnetisation 
mg calculated within linear spin wave theory (LSW) for the 3D 
Heiscnberg model Eq. (0 (upper, green dots) and the square-lattice 
J1-J2 model (lower, blue dots), as a function of J2. Remaining 
parameters for Eq. [7]are taken from experiment on CaFe2As2 [25l | . 
We use mo = 1.51, Ji x = 30meV, J\ y = l5meV , J\ z = A.5meV 
and K xy = K z = 0.12meV. The solid black line shows the sigma- 
modcl prediction Eq. (5J). The divergent correction seen in the 2D 
J1-J2 model for J2/J1 — > + 1/2 is cut off by the gap spin A a and 
3D spin-wave dispersion. 



In the highly frustrated region v y — > 0, Eq. © pro- 
vides a finite bound, 5ms < m A 3 a6c/(167r 2 xj_A cr ), on 
the maximum correction to the ordered moment from 
low-energy spin waves. To illustrate how this works, in 
Fig. [3]wc compare the predictions of the nonlinear sigma 
model, Eq. ([2]), and the Heisenberg model, Eq. (J?)), for 
the sublattice magnetization, mj, as a function of J2 
- and thereby v y . Remaining parameters for Eq. J7J) 



are taken from experiments on CaFe2As2 25j. Following 
LDA calculation^, we set the bare moment mo = 1.51^b- 
A constant offset 8ms = — 0.3/is is added to Eq. (|9|) to 
correct for high energy spin waves. The agreement is ex- 
cellent for a wide range of Ji. Even at the maximally 
frustrated point, the correction 5ms ~ 0.5/i_b is smaller 
than the 5ms ~ 0.7/zb needed to explain the discrepancy 
with experiment. We anticipate that this conclusion will 
hold for any spin model with realistic parameters^, and 
conclude that the failure of LDA to accurately describe 
the size of the ordered moment lies in high-energy elec- 
tronic correlation effects, not the zero point motion of 
low-energy spin waves. 

In conclusion, magnetic excitations in Fe pnictides are 
well-described by a simple two-fluid phenomenology in 
which gapped, three-dimensional spin waves co-exist with 
gapless, but incoherent particle-hole pairs. These two flu- 
ids can be treated as independent. This phenomenology 
is blind as to microscopic details of the real materials, 
but seems to fit naturally with spin-density wave calcu- 
lations which assign magnetism and mctallicity to differ- 
ent, weakly coupled, sheets of the Fermi surface. 
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